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1.  INTRODUCTION 

There  has  been  a  continuing  interest  in  the  use  of  the  VHF  band  to  provide 
communication  channels  which  supplement  those  of  the  more  conventional  but 
congested  HF  band.  The  availability  of  these  supplementary  channels  for  large 
range  transequatorial  circuits  was  first  reported  by  radio  amateurs  and  studies 
of  the  phenomenon  were  largely  experimental.  The  physical  mechanisms  involved 
in  producing  these  channels  have  been  and  still  are  poorly  understood.  It  is 
the  purpose  of  this  report  to  lay  down  the  elements  of  a  simple  theory  to 
explain  the  large  signal  strengths  observed  on  these  circuits  and  to  describe 
the  results  of  a  computer  model  used  to  calculate  them.. 

A  qualitative  review  of  a  variety  of  proposed  mechanisms  involved  in  producing 
these  long  distance  circuits  appeared  in  1969(ref.l).  This  review  also  contains 
an  excellent  historical  background  and  we  will  not  attempt  to  repeat  that  here. 
There  were  essentially  four  possible  ways  in  which  these  circuits  could  be  set 
up. 

(1)  Refraction  in  an  undisturbed  F  region  plasma.  The  occurrence  of  bulges 
in  the  ionization  densities  on  each  side  of  the  equator  is  assumed  to  provide 
sufficient  gradation  in  the  refractive  index  to  allow  this  to  take  place. 

(2)  Guidance  by  steep  field  aligned  gradients  in  electron  density.  The  ray 
path  is  assumed  to  be  given  a  curvature  equal  to  the  field  curvature  so  that 
guidance  occurs  along  steeply  varying  gradients  of  ionization. 

(3)  Scatter  from  F  region  irregularities.  An  ensemble  of  scatterers  at 
ionospheric  height  cannot  return  sufficient  power  to  the  opposite  hemisphere 
unless  it  is  shaped  to  provide  a  degree  of  coherence  that  is  deemed  to  be 
unlikely. 

(4)  Double  scattering  between  irregularities.  Scattering  from  irregularities 
on  one  side  of  the  equator  to  a  similar  set  in  the  opposite  hemisphere  provides 
the  correct  geometry  but  considerable  debate  exists  over  the  magnitude  of  the 
received  power. 

None  of  these  mechanisms  was  entirely  adequate  to  explain  the  phenomena  of 
transequatorial  propagation  (TEP)  and  since  the  study  was  of  a  semi-empirical 
nature  no  firm  conclusions  could  be  reached.  An  attempt  has  since  been  carried 
out  to  determine  the  physical  extent  of  field  guided  F  region  irregularities  at 
the  equatorial  peak  by  illuminating  them  with  a  sideways  looking  VHF  radar (ref. 2) 
but  the  nature  of  the  irregularities  remained  obscure. 

More  recently  large  scale  regions  of  depleted  ionospheric  plasma,  called  bubbles, 
have  been  discovered  and  observed  from  the  ground(ref . 3)  and  from  satellite 
(ref. 4).  A  description  of  TEP  of  radio  waves  at  VHF  has  also  appeared 
attributing  it  to  guidance(ref . 5)  along  field  aligned  tubes  of  depleted 
ionization  resulting  from  the  post  sunset  creation  of  equatorial  plasma  bubbles. 

A  large  ensemble  of  these  bubbles  were  then  imagined  to  contribute  the  observed 
characteristics  but  some  unusual  and  dubious  properties  were  ascribed  to  the 
waveguide  tubes. 

In  this  report  we  consider  the  properties  of  a  curved  idealized  dielectric 
discontinuity  formed  at  the  upper  edge  of  a  rising  plasma  bubble.  The  equations 
for  reflection  and  guidance  of  electromagnetic  waves  along  the  discontinuity  are 
programmed  onto  a  high  speed  digital  computer  and  profiles  for  the  received  power 
at  points  on  the  surface  of  the  earth  in  the  receiving  hemisphere  are  derived. 

A  comparison  with  some  experimental  results  for  a  Darwin  to  Yamagawa  circuit  is 
also  included. 


2.  THE  MODELS 


The  situation  is  illustrated  in  figure  1  which  shows  a  transmitter  at  point  A  in 
one  hemisphere  and  a  receiver  at  point  B  in  the  opposite  hemisphere.  The 
geographic  equator  is  located  at  G  and  the  geomagnetic  equator  is  at  position  M. 
The  magnetic  field  structure  is  assumed  to  be  identical  with  that  of  a  magnetic 
dipole  located  at  the  centre  of  the  earth  so  that  field  lines  are  symmetrical 
about  the  geomagnetic  equator. 

In  the  situation  we  have  here  a  dielectric  discontinuity  is  assumed  to  be  exactly 
aligned  with  the  particular  field  line  coincident  with  the  top  edge  of  a  rising 
plasma  bubble.  The  plasma  bubble  represents  a  region  of  depleted  ionization 
(refs. 3,4,6)  for  which  we  take  the  refractive  index  as  being  identically  one. 
Where  the  ionization  is  not  totally  depleted  then  the  refractive  index  is  only 
close  to  one  but  detailed  studies  of  this  marginal  effect  are  left  to  a  later 
report. 

Since  the  complete  modal  solution  for  this  problem  is  extremely  difficult,  if  not 
impossible,  then  it  is  necessary  to  resort  to  simplified  methods.  One  method 
that  we  can  use  in  this  situation  is  based  on  some  simple  concepts  from  geometric 
optics.  We  discuss  the  ray  approximations  in  some  detail  in  the  following 
sections. 

2.1  The  Ionospheric  Models 

It  is  shown  in  Appendix  I  that  the  field  lines  near  the  equator  assume  a  shape 
very  close  to  elliptical  (within  1%).  Furthermore  it  has  also  been  shown 
that  elliptical  multimode  dielectric  waveguiding  structures (refs . 7 ,8)  can  be 
well  represented  using  an  extended  form  of  geometric  optics  which  allows  for 
tunnelling  from  curved  dielectric  surfaces. 

An  appropriate  model  for  the  ionospheric  electron  density  is  the  Chapman  model 
described  in  Appendix  II.  The  change  in  electron  density  is  related  to  a 
change  in  refractive  index  as  given  in  Appendix  III.  Combing  these  two 
results  one  can  readily  see  that  the  refractive  index  discontinuity  varies 
with  height  above  the  earth's  surface.  The  discontinuity  gives  rise  to  a 

critical  angle  0^  for  rays  incident  from  the  more  dense  medium.  We  define 


0  =  arcos  —  (1) 

c  n:  v  J 

where  nt  is  the  refractive  index  of  the  more  dense  medium  (=1  in  this  case) 
and  n:  is  the  refractive  index  of  the  ionized  medium. 

Classical  geometric  optics  predicts  that  rays  incident  at  less  than  0 are 

totally  internally  reflected  while  rays  incident  at  greater  than  6  are 

refracted  according  to  Snell's  and  Fresnel's  Laws.  Since  the  electron 
density  varies  with  height,  then  the  critical  angle  is  also  height  dependent 
and  varies  along  the  magnetic  flux  line  which  contains  the  discontinuity. 

However,  in  this  case  because  of  the  curvature  of  the  magnetic  field  line  we 
must  extend  the  concepts  of  geometric  optics  to  include  the  phenomenon  of 
electromagnetic  tunnelling.  We  discuss  this  in  detail  in  subsequent  sections. 


At  heights  below  100  km  the  discontinuity  is  considered  to  be  so  small  as  to 
be  negligible. 


-  3  - 


ERL-0181-TR 


2.2  The  Ray  approximations 
2.2.1  Ray  concepts 

An  electromagnetic  wave  propagating  in  free  space  has  a  well  defined 
direction  of  propagation  specified  by  wave  vector  k  normal  to  the  electric 
and  magnetic  field  vectors,  E  and  H  respectively.  The  vectors  E,  H  and  k 
form  a  right-handed  triad.  In  a  uniform  homogeneous  medium,  no  aif ficulty 
exists  in  the  specification  of  the  direction  of  the  wave  and  the  concept  of 
a  ray  is  easily  managed.  A  ray  indicates  the  direction  of  the  local  plane 
wave  vector  k,  which  is  also  the  direction  of  power  propagation.  Maxwell's 
equations  are  exactly  satisfied  and  the  direction  of  the  wave  is  simply  the 
direction  of  a  ray. 

In  an  inhomogeneous  medium,  where  the  electromagnetic  fields  are  not  in 
general  specified  by  a  single  plane  wave,  the  direction  of  a  ray  is  no 
longer  so  simply  identified  and  it  is  necessary  to  specify  it  asymptotically 
as  the  short  wavelength  (X  *  0)  approximation  for  the  local  plane  wave 
direction.  It  is  used  to  indicate  the  direction  of  k  and  the  magnitude, 
given  by  the  relation 


Ik! 


k  = 


2nn 
X  * 


(2) 


is  taken  to  the  limit  k-*  ~.  X  is  the  wavelength  of  light  in  vacuum.  The 
ray  path  is  then  always  defined  by  the  eikonal  equation. 


Vn 


(3) 


where  s  is  the  distance  measured  along  the  path  and  is  the  position 
vector  for  a  point  on  the  ray  path.  n  is  the  refractive  index  of  the 
medium  at  the  point  under  consideration  and  d$?/ds  is  interpreted  physically 
as  a  unit  vector  tangent  to  the  ray  path.  This  is  the  asymptotic  theory 
known  as  geometric  optics.  The  ray  direction  defined  in  the  short  wave¬ 
length  approximation  is  specified  everywhere  including  caustics  and  points 
of  reflection. 

Indiscriminate  use  of  geometric  optics  leads  to  strange  contradictions.  It 
predicts  the  field  intensity  at  a  focus  or  caustic  to  be  infinite  which  is 
obviously  ridiculous.  The  field  distribution  near  caustics  is  not 
properly  calculated  using  geometric  optics  and  in  order  to  allow  the 
advantages  of  the  simplicity  of  geometric  optics  to  be  exploited  everywhere 
we  aim  to  extend  geometric  optics  to  accommodate  these  effects.  A  similar 
approach  has  already  been  used  by  Keller(ref . 9)  in  the  Geometric  Theory  of 
Diffraction  which  has  proved  highly  successful. 

2.2.2  Ray  tube  concepts 

Consider  the  situation  shown  in  figure  2  where  two  rays  are  drawn  through 
adjacent  points  on  a  wavefront  within  (a)  a  graded  index  and  (b)  a  step 
index  dielectric  waveguide  forming  a  characteristic  flux  tube  as  shown  in 
each  case.  The  flux  or  ray  tube  can  be  traced  as  shown  through  an  optical 
medium  by  tracing  each  of  the  rays  individually. 
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The  region  of  depleted  ionization  we  can  regard  as  being  the  core  of  a 
dielectric  waveguide  while  the  region  above  it  is  exterior  to  the  guide  and 
has  commonly  been  referred  to  as  the  cladding.  When  the  core  is  non 
absorbing,  power  in  the  ray  tube  travels  parallel  to  the  ray  direction  and 
can  only  be  lost  to  the  cladding  at  reflections  from  the  core  cladding 
interface  or  at  caustics  where  the  asymptotic  description  of  geometric 
optics  is  inadequate. 

2.2.3  A  model  for  attenuation 

The  definition  of  a  ray  allows  us  to  associate  with  a  ray  power  travelling 
in  the  same  direction.  Strictly  speaking,  the  terminology  should  refer  to 
either  a  small  ray  tube  or  flux  tube  but  it  is  simpler,  physically 
appealing  and  unambiguous  to  merely  refer  to  it  as  the  power  of  the  ray. 

In  a  guiding  structure  the  local  plane  wave  fields  of  a  family  of  rays,  all 
with  the  same  ray  invariants  and  following  similar  ray  paths,  is  the 
geometric  optics  equivalent  of  a  mode (ref . 7) .  In  order  to  develop  a  ray 
model  for  attenuation  it  is  convenient  to  introduce  the  concept  of  a  ray 
transmission  coefficient.  If  the  refractive  index  of  the  core  is  nj  and 

of  the  surrounding  medium  is  n2  with  n2  <  ni  ,  then  the  ray  suffers 
reflections  within  the  core  of  the  guide  losing  a  fraction  of  its  power  at 
each  reflection.  An  equivalent  situation  exists  for  graded  media,  where 
a  ray  loses  power  each  time  it  touches  the  outer  caustic  or  turning  point 
along  the  ray  path  in  the  core. 

The  ray  power  transmission  coefficient  T  is  defined  as 

j  -  Power  lost  from  the  incident  ray 

Power  in  the  incident  ray  ^  J 

The  value  of  T  represents  the  fraction  of  power  lost  at  each  reflection  and 
for  an  asymmetrical  guide,  such  as  we  have  here,  it  can  vary  from  one 
reflection  to  the  next.  This  definition  is  equivalent  to  the  definition 
of  the  classical  ray  transmission  coefficients  of  Fresnel. 

The  model  for  ray  attenuation  is  formulated  as  follows.  Consider  the 
situation  where  a  ray  carrying  initial  power  PQ  at  z  =  0  undergoes  a 

fractional  loss  Ti  at  the  first  reflection,  then  after  this  reflection 
there  will  be  an  amount  of  power  Pi  left  travelling  along  the  ray,  where 

P,  =  PQ(1  -  T,)  (5) 

After  a  second  reflection  we  have  power  P2  left  in  the  ray, 

P2  =  PQ(1  -  T,  )  (1  -  T2)  (6) 


where  T2  denotes  the  loss  at  the  second  reflection, 
the  remaining  power  is 


After  n  reflections 


P  =  P 
n  o 


Cl  -  T.) 
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If  the  distance  to  the  first  reflection  is  z  ,  and  the  distance  between 

I’l 

the  i  -  l'th  reflection  and  the  i ' th  reflection  is  z  .  then  the  distance 

Pi 

travelled  along  the  guide  by  the  ray  to  the  n'th  reflection  is 


n 


i=l 


It  is  convenient  to  express  the  power  variation  along  the  guide  in  the  form 
of 


Recasting  (9)  as 


P(z)  =  P(o)  exp  (  -7rz) 


-v 


=  In 


P(z) 


and  substituting  for  P(z)  and  zn  for  z  gives 


or 


'7R  )  zPi 


In (1  -  T.) 


i=l 


i=l 


Z;  ln(1  "  Ti) 


n 


i=l 


Only  when  T\  «  1  does  (12)  reduce  to  the  simple  form 


n 


i-1 


(9) 


(10J 


01) 


(12) 


(13) 
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The  waveguide  loss  between  the  transmitter  A  and  the  receiver  B  is 
therefore 


n 

ln(l  -  T.)j  (14) 

i»l 

if  we  chose  to  ignore  small  effects  such  as  the  Goos-Hanchen  shiftfref . 10) . 

This  then  becomes  a  powerful  method  for  solving  large  numbers  of  dielectric 
waveguide  problems. 

2.3  A  model  for  received  power 

So  far  we  have  examined  some  simple  two-dimensional  models  needed  for  a 
geometric  optics  based  solution  for  the  waveguiding  problem.  We  first  of  all 
need  to  extend  these  concepts  to  three  dimensions  before  we  can  calculate  the 
received  power  at  B. 

The  power  received  by  an  antenna  P  is  given  by 

K 


P(z)  =  P(o)exp 


I'E 


P  =  power  density  (D)  x  effective  area  (A) 
K 


where  the  aerial  capture  area  is 


A  = 


(IS) 


(16) 


Gr  is  the  gain  of  the  receiving  antenna  and  X  is  the  wavelength. 

The  power  density  within  a  flux  tube  is 

D  _  fraction  of  transmitted  power  in  flux  tube 

area  of  tube 

For  a  small  flux  tube  the  fraction  of  the  total  solid  angle  subtended  is 


(17) 


d0c V> 
4ir 


and  substituting  this  in  (17)  and  using  (15)  and  (16)  we  get  the  conventional 
free  space  formula 
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Where  guidance  is  involved  then  the  cross  sectional  area  of  the  tube  is  l  Rd0 
and 


Pq  d0/47r 


—  .  t.  R 
4  v 


where  1  is  the  width  of  the  flux  tube  and  d0  is  the  angular  width  of  the  tube 
in  the  dimension  in  which  guidance  occurs  as  illustrated  in  figure  3.  It  is 
important  to  note  that  this  implies  1  dependence  in  the  orthogonal  dimension, 

R 

ie  no  guidance  in  the  d <p  direction.  If  there  is  curvature  in  the  plane 
normal  to  the  xy  plane  defined  on  figure  10  then  additional  focussing  would 
be  achieved.  That  avenue  has  not  yet  been  investigated.  As  far  as  the  ray 
tracing  procedure  is  involved  we  see  that  it  is  necessary  only  to  calculate 
l.  If  guidance  also  occurred  in  the  orthogonal  direction  then  it  would  also 
be  necessary  to  calculate  the  depth  of  the  flux  tube  as  well. 

If  power  is  lost  from  the  flux  tube  by  either  absorption  or  tunnelling(ref . 11) 
then  the  power  received  must  be  attenuated  according  to  equation  (14)  since  it 
is  the  total  power  only  which  remains  constant  not  just  the  propagated  power. 
This  extension  of  geometric  optics  enables  us  to  take  account  of  curvature 
losses  in  the  same  way  as  Keller(ref . 9)  takes  account  of  diffraction  terms. 

We  therefore  have  a  method  of  calculating  the  received  power  within  each 
particular  flux  tube  enhanced  by  the  gains  of  the  receiver  and  transmitter 
antennae.  As  the  number  of  flux  tubes  increases  the  power  density  approaches 
its  limiting  value. 

If  the  transmitter  situated  at  A  has  the  antenna  polar  diagram  specified  in 
Appendix  IV  then  the  power  transmitted  into  the  ionosphere  can  be  described 
in  terms  of  ray  concepts ,  with  the  power  associated  with  each  ray  being 
determined  from  the  antenna  polar  diagram.  In  what  follows  the  gains  of  the 
antennas  will  be  normalised  to  a  standard  dipole  (G  =GD=1 . 6)  but  still 

including  the  effects  of  earth  reflection  at  both  A  and  B  according  to 
Appendix  IV. 

It  is  straightforward  to  calculate  the  points  of  intersection  of  each  ray 
with  the  magnetic  field  line  and  its  locus  in  between  until  it  intersects  the 
surface  of  the  earth  in  the  opposite  hemisphere.  The  geometry  for  the 
intersection  of  a  ray  and  the  magnetic  field  line  is  written  down  in 
Appendix  V. 

At  each  point  of  intersection  the  ray  transmission  coefficient  T  is 
calculated  using  the  Generalized  Fresnel’s  Law(ref.ll)  for  reflection  from  a 
curved  dielectric  interface.  This  law  incorporates  the  geometric  optics 
extensions  needed  to  account  for  tunnelling  transmission.  Then  summing  over 
all  flux  tubes  originating  from  the  transmitter,  contours  of  the  received 
power  for  points  near  B  can  be  derived. 
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3.  RESULTS 

The  equations  described  in  the  previous  section  were  programmed  onto  a  high-speed 
digital  computer  to  calculate  the  power  received  at  B  from  a  transmitter  at  A. 

The  power  was  normalised  with  respect  to  1  kW  of  transmitter  power  and  standard 
dipoles  at  each  end  of  a  circuit  between  Darwin  and  Yamagawa  whose  coordinates 
are  listed  in  Table  1  below. 

The  powers  received  in  dB  with  respect  to  one  milliwatt  are  plotted  in  figures 
4  to  8  against  range  from  the  geomagnetic  equator  for  wavelengths  of 
X  =  3.0,  4.0,  4.5,  5.0  and  6.0  m.  The  maximum  height  of  the  magnetic  field 

line  coincident  with  the  dielectric  discontinuity  is  taken  to  be  600  km  and  the 
critical  frequency  which  is  determined  by  the  ionization  density  is  taken  to 
correspond  with  the  predictions  of  the  Ionospheric  Prediction  Service  (IPS) 

Series  U  publications  for  March  to  April  1970.  The  position  B  corresponds  to 
the  range  of  2170  km  from  the  equator. 

Experimental  results  for  this  circuit  have  been  analysed  by  Kuriki  et  al 
(refs. 12,13)  and  plotted  against  frequency.  Their  results  were  averaged  over 
two  periods  corresponding  to  the  autumnal  and  vernal  equinoxes  for  a  period  of 
several  months  during  1970.  We  estimate  the  average  critical  frequency  for 
March  to  April  as  f£  =  10.2  Mc/s  and  for  September  to  October  as  f  =  10.8 

Mc/s.  The  results  show  an  absolute  agreement  to  within  6  dB  with  the  results 
of  our  idealised  computer  model  over  the  whole  range  of  frequencies  50  to  100 
Mc/s.  The  comparison  is  shown  plotted  in  figure  9. 

Preliminary  investigations  show  that  the  maximum  height  of  the  field  line 
required  to  give  these  results  is  about  50%  higher  than  the  height  of  the  peak 
in  the  ionization  density.  If  it  was  lower  then  the  magnitude  of  the 
dielectric  discontinuity  on  the  walls  would  be  insufficient  to  maintain  guidance 
while  if  it  was  higher  then  the  ionization  at  the  top  of  the  field  line  would  be 
insufficient. 

As  yet  the  reason  for  the  bulge  in  the  experimentally  observed  results  near 
88  MHz  has  remained  obscure  but  this  bulge  was  not  observed  in  the  results  of 
Bowen  et  al(ref. 14)  or  Nielson(ref . 1)  in  their  results.  However,  the  variation 
of  received  power  with  frequency  is  not  linear  as  proposed  by  Kuriki  et  al(ref. 13) 
even  though  over  a  limited  range  of  frequencies  it  may  appear  to  be.  He  used  a 
theory  based  on  scattering  from  ionospheric  irregularities  but  needed  to  make 
questionable  assumptions  regarding  the  nature  of  scattering  from  these 
irregularities  and  their  frequency  of  occurrence.  His  scattering  theory 
produced  the  straight  line  with  F-10  dependence  superimposed  on  figure  9. 

In  contrast  the  calculations  used  in  this  model  attribute  the  power  to  guidance 
rather  than  scattering  and  rely  merely  on  the  existence  of  field  aligned 
irregularities  that  are  likely  to  be  associated  with  rising  plasma  bubbles. 
Discrepancies  between  observed  and  calculated  results  could  easily  be  due  to 
imperfections  in  the  shape  or  continuity  of  the  elliptical  surface  assumed  in 
this  model.  Alternatively,  if  the  ionization  within  the  bubble  was  not  totally 
depleted  then  a  slightly  smaller  critical  angle  would  be  associated  with  points 
on  the  ellipse  and  this  would  reduce  the  received  power. 

Changes  in  the  critical  frequency  of  6%  observed  between  the  autumnal  and  vernal 
equinoxes  produces  only  a  dB  change  in  the  received  power  for  the  Darwin/ 
Yamagawa  circuit  at  X  =  6m  using  this  geometry  rather  than  the  3  dB  observed 
experimentally. 


I 
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Since  the  results  of  this  initial  investigation  are  so  close  to  the  experimental 
results  further  analysis  appears  to  be  warranted.  In  a  subsequent  study  we  will 
perturb  all  the  parameters  associated  with  this  model  to  determine  their  allowable 
range  and  to  seek  an  experimental  method  of  confirming  their  behaviour. 


4.  CONCLUSIONS 

A  simple  model  for  the  received  power  on  a  transequatorial  circuit  described  by 
figure  1  has  been  proposed  which  may  explain  night-time  radio  wave  propagation  in 
the  ionosphere.  The  mechanism  has  been  ascribed  to  guidance  along  field  aligned 
irregularities  associated  with  rising  plasma  bubbles  in  the  ionosphere.  Good 
agreement  has  been  demonstrated  between  experimental  results  obtained  on  a  circuit 
between  Darwin  and  Yamagawa  and  theoretical  results  for  an  ideally  formed  guiding 
irregularity. 

Further  investigation  is  required  to  optimise  the  model,  to  determine  an 
experimental  technique  of  validating  it  and  to  evaluate  the  effect  of  imperfec¬ 
tions  in  the  guiding  surface. 


ERL-0181-TR 


10  - 


REFERENCES 


No.  Author 


Title 


1 


2 


3 


4 


5 


6 


7 


8 


9 

10 


11 


Nielson,  D.L. 


Booker ,  H . G .  and 
Ferguson,  J.A. 


Woodman,  R.F.  and 
La  Hoz,  C. 


Dyson,  P.L.  and 
Benson,  R.F. 


Heron,  M.L.  and 
McNamara ,  L . F . 


McClure,  J.P.  et  al 


Love,  J.D.  et  al 


Winkler,  C. 

Keller,  J.B. 

Love,  J.D.  and 
Winkler,  C. 

Snyder,  A.W.  and 
Love,  J.D. 


"Long  range  VHF  propagation  across  the 
geomagnetic  equator". 

Stanford  Res.  Ins.  Report,  March  1969 

"A  theoretical  model  for  equatorial 
ionospheric  spread  F  echoes  in  the  HF  and  VHF 
bands". 

J.  Atmos.  Terr.  Phys.,  Vol.  40,  p  803, 

January  1978 

"Radar  observations  of  F  region  equatorial 
irregularities" . 

J.  Geophys.  Res.,  Vol.  81,  p  5447, 

November  1976 

"Topside  sounder  observations  of  equatorial 
bubbles". 

J.  Geophys.  Res.,  Vol.  5,  p  795, 

September  1978 

"Transequatorial  VHF  propagation  through 
equatorial  plasma  bubbles". 

Radio  Sci.  Vol.  14,  p  897,  Sept-Oct  1979 

"Plasma  bubbles  and  irregularities  in  the 
equatorial  ionosphere". 

J.  Geophys.  Res.,  Vol.  82,  p  2650,  July  1977 

"Rays  and  modes  in  step  index  multimode 
elliptical  fibres". 

Micro.  Opt.  Accous.,  Vol.  3,  p  231, 

November  1979 

"Theory  and  application  of  ray  transmission 
coefficients  in  multimode  optical  fibres". 
Ph.D.  thesis,  Aust.  Nat.  Univ.  1978 

"Geometric  theory  of  diffraction". 

J.  Opt.  Soc.  Am.,  Vol.  52,  p  116,  July  1962 

"Goos-Hanchen  shift  for  leaky  rays  on  step 
index  waveguides". 

Electron.  Letts.,  Vol.  14,  p  379,  June  1978 

"Reflection  at  a  curved  dielectric  interface 
-  electromagnetic  tunnelling". 

IEEE  Trans.  WTT-23,  p  134,  January  1975 

"Propagational  mode  deduced  from  signal 
strengths  in  the  VHF  band  on  the  trans¬ 
equatorial  path". 

J.  Radio  Res.  Labs.,  Vol  19,  p  175,  June  1972 


12 


Kuriki,  I.  et  al 


11 


I-RL-0181-TR 


No. 

Author 

Title 

13 

Kuriki,  I.  et  al 

"Some  characteristics  of  VHF  propagation 
between  Australia  Japan". 

5th  International  Symposium  of  Equatorial 
Aeronomy,  Townsville,  August  1976 

14 

Bowen,  E.D.  et  al 

"VHF  characteristics  of  the  transequatorial 
ionosphere". 

J.  Geophys.  Res.,  Vol.  73,  p  2469,  April  1978 

15 

Winch,  R.P. 

Electricity  and  Magnetism,  Prentice  Hall,  1959, 
p  350 

16 

Budden ,  K .  G . 

Radio  Waves  in  the  ionosphere,  Cambridge 
University  Press,  1961 

ERL-0181-TR 


12 


k  ,  k  .  k  ,  k 
yi  yz  ya  y4 


Ill  ,  112 


NOTATION 


intercept  of  ray  with  y  axis 
variable  defined  in  Appendix  V 
critical  frequency 
plasma  frequency 

maximum  height  of  field  line  chosen 


wavenumber 


component  of  wave  vector  I  kl 
ray  slope 

incident  ray  slope 

slope  of  ellipse  at  point  of  reflection 
reflected  ray  slope 
refractive  indices 

polar  coordinate  for  point  on  field  line 

minor  semi  axis  of  ellipse 

major  serai  axis  of  ellipse 

ionospheric  scale  height 

antenna  height 

electron  density 


radius  of  earth 


ray  transmission  coefficient 


Fresnel  coefficient 


coefficient  of  °c  Chapman  layer 
dielectric  permittivity 
wavelength  of  source 
critical  angle 
angle  of  incidence 

radius  of  curvature  of  ellipse  at  point  of 
reflection 
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phase  path  difference 

ray  direction  from  transmitter 

variable  defined  in  Appendix  IV 
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APPENDIX  I 

STRUCTURE  OF  THE  EQUATORIAL  FIELD 

The  equation  of  a  metnetic  field  line  from  a  point  magnetic  dipole  located  at  the 
centre  of  the  earth  is  given  by  the  equation (ref . 15) 

T 

2a  =  constant  (1.1) 

COS  P 


where  r,  0  are  the  polar  coordinates  of  a  point  on  the  field  line  as  shown  in 
figure  10.  The  intersection  of  a  field  line,  which  passes  through  a  maximum 
height  h  at  the  geomagnetic  equator  (G) ,  with  the  surface  of  the  earth  (radius 
R)  occurs  where  latitude  0  is  given  by 


6  =  ±arccos  J  (I>2) 

In  the  local  plane  cartesian  coordinates  Oxy  position  of  the  intersection  is  at 
(Rsin0  ,  Rcosfl  -  R/2). 

If  the  discontinuity  is  assumed  to  begin  at  conjugate  points  on  opposite  sides 
of  the  equator  S,  S'  at  a  height  of  100  km  then  in  local  cartesians  Oxy  the 
points  S,  S'  using  equation  (1.2)  are  at  (R  +  100)  sin0 ' ,  (R  +  100)  cos0 ' ,  - 
R/2,  where  0'  is  the  latitude  of  intersection  of  the  field  line  contour  and  the 
contour  100  km  above  the  earth's  surface.  Fitting  an  ellipse  to  the  points  G, 
S,  S'  gives  a  best  fit  curve  having  the  equation 


_ X _  _■) 

(R  +  100)  sin0^. 

^  j 

J  i  .  /(R  +  1OO)cos0  '  -  R/2 

»  \ - (TTh)72 - / 

♦  |  _ L _ 

V  (R  +  h)/2 


1 

(1.3) 


In  the  region  of  the  equator  the  error  involved  in  assuming  the  equation  of  the 
field  line  to  be  equation  (1.3)  for  h  =  500  km  is  less  than  1%  compared  to 

equation  (I . 1)  over  the  whole  of  the  curve  above  the  ground. 
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APPENDIX  II 

MODEL  FOR  ELECTRON  DENSITY 

An  appropriate  electron  density  model  is  the  Chapman  model  defined  by  = 
Reference  16.  The  frequency  f  is  defined  in  terms  of  the  critical  frequency 
fc  by 

fN  =  fc  6XP  (<X(1  "  Z  '  6"Z)) 


where 


z 


h  -  h  max 
H 


a  =  0.5  for  an  a  Chapman  layer 


h  max  is  the  height  of  the  electron  density  peak. 

A  plot  of  height  against  ionization  density  for  varying  values  of  H  with 
f  =  12  MHz,  h  max  =  400  km  appears  as  Figure  11.  We  choose  H  =  60  km  to 

give  the  distribution  of  Booker (ref . 2) . 

The  electron  number  density  N  is  then  given  by  the  formula 

N  =  1.24  x  10"8  f?t 

N 


J 


APPENDIX  III 


CHANGE  IN  REFRACTIVE  INDEX  AT  AN  IRREGULARITY 


The  fractional  change  Ac  in  the  dielectric  permittivity  e  in  an  ionized  medium  is 
related  to  the  change  in  ionization  density  An  by  the  relation 


_1_  2 

*2 


X4 


(AN)2 


where  r  is  the  classical  radius  of  the  electron 
e 


(2.82  x  10_,s  m) 


Strictly,  X  should  be  the  wavelength  taking  into  account  the  ionization  density 
N,  but  for  VHF  propagation  in  this  model  we  will  take  X  to  be  the  free  space 
wavelength.  At  the  irregularity  therefore  we  will  see  a  change  in  the  magnitude 
of  the  refractive  index  discontinuity  with  height  and  a  corresponding  change  in 
the  critical  angle  0  where 


6  =  arcos  — 

c  n2 


where  nj ,  m  are  the  refractive  indices  on  either  side  of  the  discontinuity. 

In  the  case  of  a  totally  depleted  plasma  bubble  rising  in  the  ionosphere  we  put 

m  =1. 

If  we  arrange  =  1  then  n2  is  determined  and  hence  0  . 

N  C 


17 


ERL-0181-TR 


APPENDIX  IV 

ANTENNA  POLAR  DIAGRAMS 

The  polar  patterns  characterising  either  a  tiansmitting  or  receiving  antenna  at 
height  Hi  above  the  ground  are  assumed  to  be  identifical .  The  received  field 
is  given  by  the  sum  of  the  direct  and  reflected  wave  amplitudes  as  depicted  in 
Figure  12.  At  low  angles  we  can  assume  that  the  reflection  coefficient 
of  the  surrounding  earth  is  -1  and  if  the  wave  amplitude  is  A  then  the  total 
amplitude  received  assuming  unity  antenna  gain  is 


=  Asinwt  -  A  sin(wt  +  <t> ) 


(IV. 1) 


where 


and 


This  reduces  to 


X  =  operating  wavelength 


A  =  angle  of  propagation. 


A|sinwt(l  -  cos$)  -  coswt  sin^j 
=  A'  sin(wt  -  4*) 


(IV. 2) 


where 


<t>  = 


(IV. 3) 


The  magnitude  of  the  power  from  a  transmitting  antenna  in  direction  A  is  then 


A'2  = 


A2  (2  -  2cos 0) 


(IV. 4) 
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APPENDIX  V 


THE  GEOMETRY  OF  REFLECTION  FROM  AN  ELLIPTICAL  SURFACE 

Simple  geometry  as  shown  in  figure  13  is  used  to  determine  the  angles  of 
incidence  and  reflection  at  a  dielectric  discontinuity.  If  the  equation  of  the 
ellipse  is 


A2 


O' .  -  YI2 

bt 


i 


and  of  the  incident  ray  is  y  =  mx  +  b  where  b  is  determined  by  the  initial 
conditions  then  the  point  of  intersection  is  given  by  the  solution 

-A2  me  +  ABv'B2  -  c2  +  A2m2 

x  =  - (FT??) 


where 


c  =  b  -  Y 


The  slope  of  the  ellipse  at  the  point  of  intersection  is  m2  given  by 


-xB 

1,12  =  T2 

If  the  angle  of  incidence  measured  with  respect  to  the  reflecting  surface  is 
0 i  *  arctan  m2  -  arctan  mi  and  the  slope  of  the  reflected  ray  is  m3  then 

m3  =  tan  0 

r 

and 


m3  -  2  m2  ~  mi 


* 

k 
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TABLE  1.  COORDINATES  OF  DARWIN  AND  YAMAGAWA 


GEOGRAPHIC  COORDINATES 

LATITUDE  RELATIVE  TO 

GEOMAGNETIC  EQUATOR 

DARWIN 

130°  50' 

E,  12°  20'  S 

19°  32'  S 

YAMAGAWA 

130°  37' 

06"  E,  30°  12'  05"  N 

23°  O'  N 

on  of  flux  ti 
m  a  wavefron' 
ectric  wavegi 


UOO  1600  1800  2000  2200  2400  2600  2800 

eANGE (  km) 

Figure  6.  Plot  of  normalised  received  signal  strength  against  distance 
from  the  equator  for  X  =  4.5  metres  and  h  =  600  kilometres 
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Figure  7.  Plot  of  normalised  received  signal  strength  against  distance 

from  the  equator  for  X  =  5.0  metres  and  h  =  600  kilometres 


km) 

Figure  8.  Plot  of  normalised  received  signal  strength  against  distance 
from  the  equator  for  X  =  6.0  metres  and  h  =  600  kilometres 
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Figure  9 


Figure  9 


Plot  of  normalised  received  signal  strength  against  frequency 
for  both  theoretical  and  experimental  results  as  discussed  in 
the  text 
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Figure  11.  Plot  showing  relat: 
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